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ABSTRACT 



We present high spatial resolution ^^CO (J = 1 — 0) images taken by the 
Nobeyama 45m telescope toward a 48' x 48' area including the L1641-N cluster. 
The effective spatial resolution of the maps is 21", corresponding to 0.04 pc 
at a distance of 400 pc. A recent 1.1 mm dust continuum map reveals that 
the dense gas is concentrated in several thin filaments. We find that a few dust 
filaments are located at the parts where ^^CO {J — 1 — 0) emission drops sharply. 
Furthermore, the filaments have two-components with different velocities. The 
velocity difference between the two-components is about 3 km s~^, corresponding 
to a Mach number of 10, significantly larger than the local turbulent velocity in 
the cloud. These facts imply that the coUision of the two components (hereafter, 
the cloud-cloud collision) possibly contributed to the formation of these filaments. 
Since the two components appear to overlap toward the filaments on the plane 
of the sky, the collision may have occurred almost along the line of sight. Star 
formation in the L1641-N cluster was probably triggered by such a collision. 
We also find several parsec-scale CO shells whose centers are close to either the 
L1641-N cluster or V 380 Ori cluster. We propose that these shells were created 
by multiple winds and/or outflows from cluster YSOs, i.e., "protocluster winds." 
One exceptional dust filament located at the western cloud edge lies along a shell; 
it is presumably a part of the expanding shell. Both the cloud-cloud collision and 
protocluster winds are likely to influence the cloud structure and kinematics in 
this region. 

Subject headings: ISM:structure — ISM: clouds — ISM: kinematics and dynamics — 
stars: formation 
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Introduction 



Most stars form in giant molecular clouds (GMCs). In GMCs, various environmental 
effects such as large-scale flows, supernovae, and stellar feedback from young stars 
(winds, radiation, and outflows) often shape the cloud structure a nd dynamics, trigge ring 



McKee fc Ostrikei 



and suppressing the formation of the next-generation stars (e.g., 
20071 ). Recent numerical simulations of star formation have demonstrated that 
these effect s significantly influence, o r even control, the cloud evolution and star 



formation (Mac Low fc Klessen 



2004 



Krumholz et al. 



2OO9I 



2011 



Li &: Nakamura 



2006 



Nakamura &: Li 



2007 



Baner 



Vazquez- Semadeni et al 



ee et al. 



Carroll et al 



2009 



Gritschneder et al 



l2009l : 



2OIOI ). Fi ngerprints of these environmental effects have been founc 



in some star-forming 



2008 



Sandell fc Knee 



2001 



^ (e.g.. 


Bally 


1989 


: Hever et al. 


1992; 


Ballv 


2008; 


O'Del 


et al. 


Shimajiri et al. 


2008; 


Arce et al. 


2OI0I; 


Shimajiri et al. 


2011) 



However, the roles of the environmental effects in star formation remain poorly understood 
observationally. This is partly because wide-field, high spatial and/or spectral resolution 
observations, which are needed to resolve the cloud structure and kinematics in detail, are 



still limited. In parti cular, stellar 



to parsec-scale (e.g.. 



Heyer et al. 



feedbacks such as winds or outflows are often extended 



19921 ). Wide-field observations of the cloud structure 
and kinematics are needed to unveil such environmental effects. At the same time, it is 
necessary to resolve the cloud structure at a scale of "dense cores", which are the basic 
units of individual star formation 0.1 pc ~ one arcmin at a distance of 400 pc), in order 
to uncover a link between individual star formation and the environmental effects. 

In this paper, to understand how the environmental effects influence the internal 
structure and kinematics in star-forming molecular clouds, we present the results of 
wide-field ^^CO (J = 1 — 0) mapping observations toward the L1641-N region, a nearby 
active star-forming region in the Orion A giant molecular cloud complex, using the 
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Nobeyama 45 m radio telescope. Our data have high angular 21") resolution, allowing 



us to resolve spatia 



Menten et al 



2007 



structures at a scale of 0.04 pc at a dist ance of 400 pc (see e.g. 



Sandstrom et al. 



2007 



Hirota et al. 



19921 ). By comparing the ^^CO 



( J = 1 — 0) map with the 1.1 mm dust continuum map taken by the AzTEC camera on the 
ASTE telescope, we investigate how the dense gas is distributed in the parent molecular 
cloud. 



The L164 



flAUen &: Davis 



-N re gion is one of the well-studied, nearby star-forming regions 

20081 . and references therein). It lies just south the clouds known as OMC-4 



and OMC-5, making up a northern part of the L1641 molecular cloud, one of the several 
~ IO^Mq molecular clouds contained within the Orion A giant molecular cloud complex. 
The L1641 cloud is a filamentary cloud extending more than 2 — 3 degrees (15 ~ 20 pc) 
and having numerous clumpy elongated condensations with typical masses of a few tens 
to ~ IO^Mq and sizes of a few tenths to ~ 1 pc. This cloud is not forming massive stars: 
its most massive member is the B4 V Herbig Ae/Be star, HD 38023, at the position of 
(R.A., Dec.) = (5:42:21, -8:08), located at the L1641-S region. Although it has no rich 
clusters comparable to the Orion Nebula Cluster, the L1641-N region consists of several 
active star-forming regions. Optical, X-ray, and infrared surveys have revealed that most o: 



the protostars are clustered in small clusters or groups with a few xlO YSOs (IStrom et al. 



1993 



Carpenter et al 



approximately 80 YSOs ( ICalfalk fc Olofsson 



200lh. A prominent examp 



2008 



e is the L1641-N cluster which contains 



Fang et al. 



20091 ) . having the stellar 



surface density of ~ 200 pc ^, an order of magnitude larger than the average stellar surface 
density of the distributed population of YSOs. The most luminous object in t his cluster is 



IRAS 05338-0624, around 



1986 



Stanke fc Williams 



which powerful molecular outflows were discovered (IFukui et al. 



20071 ). 



Another example is the V 380 Ori cluster. V 380 Ori itself is a binary of Herbig Ae/Be 



- 5 - 



stars, illuminating the reflection nebula NGC 



(H-H) objects were discovered in this region (j Allen fc Davis 



999 flBaines et all 120061) . Many Harbig-Haro 



20081 . and references therein). 



The most fa mous and spectacula r outflows are H-H 1/2, which are associated with gigantic 



bow shocks. 



Heyer et al. 



I1992I ) discovered a parsec-scale expanding shell around the V 
380 Ori star from the ^^CO (J = 1 — 0) observations, suggesting that star formation activity 
indeed shapes cloud structure and dynamics significantly. These two clusters do not contain 
massive stars that would emit strong UV radiation and control the evolution of the region. 
The L1641-N region is expected to provide us with clues to understanding the impact of 
the current and previous star formation activity on the cloud structure and kinematics. 

The rest of the paper is organized as follows. The details of our observations and data 
are described in Sections |2] and [31 We present results of our ^^CO ( J = 1 — 0) observations 
in Section HJ and discuss the cloud structure and star formation activity of this region in 
Section [51 Finally, we summarize our conclusion in Section [61 



2. Observations 

The ^^CO {J = 1 — 0; 115.271204 GHz) observations were carried out with the 
25-element focal plane receiver BEARS on the Nobeyama Radio Observatory (NRO) 45 m 
telescope. They cover a 48' x 48' area including the L1641-N cluster, the northern part 
of the L1641 molecular cloud, from 2009 December to 2010 January. At 115 GHz, the 
telescope has a FWHM beam size of 15" and a main beam efficiency, 77, of 0.32. At the 
back end, we used 25 sets of 1024 channel auto-correlators (AGs) which have bandwidths 
of 32 MHz and frequency resolutions of 37.8 kHz. The frequency resolution corresponds to 
a velocity resolution of 0.1 km s~^ at 115 GHz. During the observations, the system noise 
temperatures were in the range between 300 to 600 K in DSB at the observed elevations. 
The standard chopper wheel method was used to convert the output signal into the antenna 
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temperatures (T^), corrected for the atmospheric attenuation. Our mapping observations 
were made by the On-The-Fly (OTF) mapping technique. We adopted a spheroidal function 
as a gridding convolution function (GCF) to calculate the intensity at each grid point of 
the final cube data with a spatial grid size of 7.5" and a velocity resolution of 0.5 km s~^. 
The final effective resolution of the map is 21", corresponding to 0.04 pc at the distance to 
Orion A of 400 pc. The rms noise level of the final map is 0.39 K in T^. 



3. Other Data 



In the next section, we compare our ^^CO {J = 1 — 0) data with 1.1 mm continuum, 
i3co (J = 1 - 0), and H^3C0+ ( J = 1 - 0) data. Here, we briefly describe these data. 

The 1.1 mm continuum data were taken toward a 1.7° x 2.3° region in the northern 
part of the Orion A giant molecular cloud comp lex with the AzTEC camera mounted on 



the ASTE 10 m telescope (jShimajiri et al. 



1201 ll ). The observations were carried out in the 
period from October to December 2008. The noise level was about 9 mJy beam^^ and the 
effective beam size was 40", about twice that of the ^^CO map, after the FRUIT imaging 
which is an iterative mappin g method to recover th e spatially-extended component. The 



detail of the data is given in 



Shimajiri et al. 



(120111 ) 



The "CO ( J = 1 - 0) data were taken by iBallv fc Langerl (119871 ) with the 7 m 



telescope of the AT&T Bell L aboratories. The orig i nal data co vers the whole Orion A giant 



molecular cloud complex (see 



Bally fc Langer 



1987 



Bally 



19891 ). much larger area than that 



of our ^^CO ( J = 1 — 0) data. The observations were carried out in the period from 1984 to 
1986. The grid size of the data was 60", about three times coarser than that of our ^^CO 
(J = 1 — 0) data. The rms noise level was about 0.3 K in T4 at the velocity resolution of 
about 0.27 km s~^. The main beam efficiency of the 7 m telescope was 0.9. 
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The H^^CO"*" ( J = 1 — 0) data were taken from the Nobeyama 45 m archival data. The 
data were obtained from December 1999 to April 2004 in the position-switching mode, with 
a grid size of 21". The rms noise level was 0.1 K in at the velocity resolution of about 



0.13 km s ^ with a main-beam efficiency of 0.51. See http://www.nro.nao.ac.jp/ for more 
detail. 



4. Results 

4.1. Global i^CO ( J = 1 - 0) and ^^CO ( J = 1 - 0) Distributions 

In Figure [H we present a ^^CO ( J = 1 — 0) total integrated intensity map toward 
the 48' X 48' area including the L1641-N cluster. The ^^CO (J = 1 — 0) emission tends 
to be optically thick in the entire molecular cloud, thus representing the spatial extent of 
the global molecular gas distribution, instead of the density distribution. For comparison, 
in Figure [T](b), the 1.1 mm continuum map taken by the AzTEC camera mounted on the 
ASTE 10 m telescope is overlaid on the ^^CO (J = 1 — 0) total integrated intensity map 
with the contours. 

The ^^CO {J = 1 — 0) integrated intensity map indicates that the parent filamentary 
cloud is roughly running from north to south. The filament appears to bifurcate in the 
southern part of the image. At the bifurcation point, the L1641-N cluster resides. The 
most intense ^^CO emission is associated with the L1641-N cluster at the position of (R.A., 
Dec.) = (5:36:18, —6:21:48). The position of the ^^CO integrated intensity peak almost 



coincides with that o: 



Stanke fc Williams 



' the bright, compact continuum source, L1641-N MMl identified by 



( 120071 ) with the SMA interferometer, although it deviates about 22" 
toward south. This is probably because the single-dish continuum peak traces the whole 
clumpy structures identified by SMA, including L1641-N MMl, as shown in Figure 6 of 
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Stanke &: Williams I (120071 ). Comparison between the 1.1 mm continuum image and the 



^^CO integrated intensity map indicates that several filamentary dense clumps are located 
near the parts where the ^^CO emission drops sharply. These parts are in good agreement 
with the edges traced by ^^CO (J = 1 — 0) emission (see Fig. [2]) and therefore, hereafter, 
we refer to these parts as the cloud edges, for simplicity. For comparison, several main 
filamentary dense clumps are labeled with alphabet A through E in Figure [T]|^b), indicating 
that the filaments A and E are located at the eastern and western cloud edges, respectively. 

In Figure we also present the ^^CO ( J = 1 — 0) integrated intensity map toward the 
same area shown in Figure [TJ Since the fractional abundance of ^^CO is smaller than that 
of ^^CO by a factor of about 40, the ^^CO emission tends to be optically thinner than the 
^^CO emission, and thus its spatial distribution is expected to represent the distribution of 
molecular gas with intermediate densities of about 10'^ cm~'^. Although the spatial extent 
of the ^^CO integrated intensity is well covered by that of ^^CO, its spatial distribution 
appears different from that of ^^CO. The ^'^CO emission tends to be concentrated in 
elongated or filamentary structures, which reasonably trace the parsec-scale filamentary 
clumps identified by the 1.1 mm dust continuum emission (see Figure |2]). The 1.1 mm dust 
continuum emission and ^^CO ( J = 1 — 0) integrated intensity also take their maxima 
nearly at the same position, i.e., toward the L1641-N cluster center, suggesting that the 
dense gas is concentrated into the L1641-N cluster-forming filamentary clump (the filament 
C). 

In Figure [31 we present a ^^CO (J = 1 — 0) peak intensity map, which indicates that 
the ^^CO (J = 1 — 0) emission drops sharply both at the eastern and western sides of the 
parent filamentary cloud. At the same places, the ^^CO and ^^CO integrated intensity 
drops sharply. A prominent feature is the existence of many arcs or elongated structures, 
most of which tend to be across the filamentary structures traced by the ^^CO ( J = 1 — 0) 
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integrated intensity map. The peak intensity map of the optical thick ^^CO (J = 1 — 0) 
emission typically represents the distribution of the excitation temperature, instead of the 
local density distribution. However, the structures traced by the ^^CO {J = 1 — 0) peak 
intensity map are sometimes seen in the ^^CO ( J = 1 — 0) integrated intensity map and is 
therefore likely to reflect the cloud density distribution in some parts. For example, some 
arcs located at the southwest part near the cloud edge are recognized in both the ^^CO 
(J = 1 — 0) peak intensity and ^'^CO integrated intensity maps (e.g., (R.A., Dec.) ~ (5:35:0, 
-6:40:0) and (5:35:30,-6:35:30)). 



4.2. Dust Filaments 

Here, we present several aspects of the parsec-scale filamentary dust clumps, by 
comparing the spatial distributions of several molecular emission lines such as ^^CO 
( J = 1 — 0) and ^^CO ( J = 1 — 0) with the 1.1 mm continuum emission. 

The dust filaments are likely to have high densities ~ 10^~^ cm~^ because the CS 
( J = 1 — 0) emission is detected toward all the dust filaments (see Figure 5(c) of Tatematsu 
et al. 1993). Th e H^^CO"*" (J = 1 — 0) emission is also detected toward the dense parts 



of the filaments fllkeda et al 



20071 ). For comparison, we show in Figure H] the H^^CO"'" 
(J = 1 — 0) integrated intensity contours overlaid on the 1.1 mm image, in which the dust 
filaments are labeled with A through E. Figure H] indicates that the strongest H^'^CO^ 
( J = 1 — 0) emission comes from the most massive filament C associated with the L1641-N 
cluster. The H^^CO^ emission peaks in this filament coincide reasonably well with those of 
the 1.1 mm continuum emission. The line-of-sight velocities of the H^^CO^ (J = 1 — 0) 
peaks are around 7 km s"^. Interestingly, the CS (J = 1 — 0) emission is strongest at Vlsr ~ 
10 km s~^, but weaker at Vlsr ~ 7 km s^^ (see Figure 5(c) of Tatematsu et al. 1993), 
although the H^'^CO^ emission is strong at around 7 km s^^. In addition, the CS emission 
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associated with the massive filament C is distributed in the wide velocity range of 6 km 
s^^ to 12 km s^^. As mentioned below, the CS emission as well as ^^CO indicates that the 
most massive filament C has two different velocity components, which may give us a clue to 
understand the formation mechanism of the filament. We note that the filament A is out of 
the area observed by 



Tatematsu et al. 



IkedaetaL 



([199J Thus, the filament contains the dense gas with densities of 10^ ^ 



(120071), but it is detected in CS (J = 1 - 0) obtained by 



cm ^. 

The ^^CO ( J = 1 — 0) velocity channel maps presented in Figure indicate that the 
velocity structure in this region can be divided into three components: blueshifted (3—6 
km s~^), main (7—9 km s~^), and redshifted ones (10—14 km s~^). For comparison, we 
present in Figure [6] a two-color image with redshifted ^^CO integrated intensity in red and 
blueshifted ^^CO integrated intensity in blue. The 1.1 mm continuum map is overlaid on 
the image with white contours. The integration ranges are 9.5 ~ 14.5 km s~^ and 3.5 
~ 6.5 km s~^ for the redshifted and blueshifted components, respectively. The two-color 
image indicates that the redshifted component is dominant in the upper half of the parent 
filament, whereas the blueshifted component is dominant in the lower half. 

The strong blueshifted component with a velocity range between 4 and 6 km s^^ 
appears in the southern part where the filament C and the L1641-N cluster reside. The 
redshifted component with a velocity of 10 to 12 km s^^ is also associated with the massive 
filament C. The presence of a couple of velocity components toward the filament C is more 
clearly seen in the position- velocity map presented in Figure [7|(c), showing that the filament 
consists of two separate components with different velocities toward the eastern part of the 
filament: one is about 6 km s~^ and the other is 9 km s~^. Th ese two-compone n ts wit h 



different velocities are in good agreement with the CS map by .Tatematsu et al. 



(119931 ) 



The filament is surrounded by the blueshifted elongated component that apparently have 
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a head-tail shape (see also Figure [6]). Another faint velocity component is seen toward the 
western part of the filament at the velocity of about 11 km s^^. This component is more 
evident in the position-velocity diagrams of the ^^CO {J = 1—0) emission (see Figure [HI^c)). 

Other dust filaments also have two-component velocity structures. The position-velocity 
diagram (Figure [7](b)) indicates that the filament A, located at the eastern cloud edge, has 
two-components with different velocities: a diffuse component is at 6 km s~^ and stronger 
emission component is at 9 km s~^. The similar velocity components are associated with 
the filaments B and D (see also Figure [5]). The two-color image indicates that the eastern 
area outside the parent cloud i s filled wit h the diffuse componen t with about 6 km s~^. 



According to 



Sakamoto et al. 



fll997h and 



Shimajiri et al. 



(1201 ll ). this diffuse component 



is distributed in the much larger area in the eastern area of the Orion A giant molecular 
cloud complex. The filament E has the two-components at 7 km s~^ and 9.5 km s (see 
Figure [Tll^e)). These velocity structures of the filaments are in good agreement with the 
position-velocity map of the ^^CO ( J = 1 — 0) emission (see Figure [8](e)). It is worth noting 
that the two filaments labeled with A and E are located near the cloud edges where the 
^^CO ( J = 1 — 0) emission drops sharply. 



4.3. Shells 

Besides the filamentary structures, shell-like structures are prominent in our ^^CO map, 
particularly in the peak intensity map. In this subsection, we describe some characteristics 
of these shell-like structures in detail. 

In the ^^CO ( J = 1 — 0) integrated intensity map (Figure [1]), a shell-like structure can 
be recognized at the position of (R.A., Dec.) ~ (5:36:9, —6:4:20). The radius of the shell is 
about 10', corresponding to about 1 pc at a distance of 400 pc. The shell is more clearly 
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seen in the CO (J = 1 — 0) integra" 



Heyer et al. 



;ed intensity map presented in Figure [2l This shell 



(I1992I ) who carried out a wide-field mapping observation in 



was first found by 

^■^CO ( J = 1 — 0) toward the Orion A giant molecular cloud complex with the FCRAO 14 m 
telescope. They found that the shell has two velocity components with a velocity difference 
of about 2—3 km s~^, implying that it has an expanding motion. They also found a number 
of large holes surrounded by expanding shells. The typical radius of the shells is 10' — 26', 
corresponding to 1—3 pc at a distance of 400 pc. They interpreted that these structures 
stem from events associated with the energetic star formation activity within the cloud. 

Similar shell-like structures can be seen in our ^^CO maps. For example, in the central 
part of the ^^CO peak intensity map, two prominent thin shell structures are seen (Figure 
[3]). We labeled these shells by Al and A2 on the peak intensity map presented in Figure [91 
There is also a shell labeled by A3, which is almost parallel to the Al and A2 shells. These 
three shells are spatially well-ordered and appear to have a common center which is very 
close to the position of the L1641-N cluster. These shells have very small thickness of about 
30". Such a thin shell may be difficult to be found in the ^'^CO map because of the spatial 
resolution. 

Another prominent shell-like structure can be recognized in the southern part of the 
CO integrated intensity maps (Figures [U and [21), where the CO emission appears weak. 
Here, we labeled three remarkable shell-like structures or arcs by Bl, B2, and B3, although 
there are several similar arcs seen in the image (Figure [9}). These arcs appear to be spatially 
well-ordered an d homocentric. Some o f them are detected in the CS (J = 1 — 0) line (see 



Figure 5(c) of iTatematsu et al.lll993l ). indicating that the dense gas is associated with 
them. We note that the common center of these shells is very close to the position of the V 
380 Ori cluster, which is just outside our ^^CO map, located at (R.A., Dec.) = (5:36:25.43, 
-6:42:57.7). 
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The shells are also clearly seen in a volume rendering image of the ^^CO ( J = 1 — 0) 
antenna temperature data presented in Figure [101 The volume rendering technique is 
suitable to find and visualize the coherent structures in the 3D space. The green color 
typically represents the areas with relatively high of about 10—15 K. We designated 
the shells Al, A2, A3, Bl, B2, and B3 in the image. These shells are coherent in the 
position-position-velocity cube, suggesting that they were created by some dynamical 
events. Several arcs, parts of the shells, can also be recognized in the ^^CO map presented 
in Figure |2j These coherent structures are remarkable since the surrounding gas is highly 
turbulent and such coherent structures would have been destroyed in a crossing time of 
the shells (about 10^ yr for the radius of 10" — 30" and velocity of 3 km s~^). Their 
existence suggests that the events producing these coherent structures play a dominant role 
in determining the cloud internal structure. In Section E3] we will discuss these structures 
in detail. 



4.4. Molecular Outflows in the L1641N cluster 



Recent numerical simulations of cluster formation have suggested that protostellar 
outflow feedback plays an important role in regulating star formation in cluster-forming 
clumps because they inject substantial amount of kinetic energy into the surround ings 



( iLi fc Nakamura 



2006 



Nakamura &: Li 



2007 



Carroll et al. 



2009 



Wang et al. 



20101 ). Here, 



we attempt to identify molecular outflows associated with the L1641-N cluster to assess the 
roles of the outflow feedback in the dynamical evolution of the cloud. The ^^CO emission 



has been successfully adopt ed to identify high-velocity components driven 



protostellar outflows, (e.g., 



2010 



Nakamura et al. 



Stanke fc Williams 



2011a 



2007 



Takahashi et al. 



) y powerful 



2008 



Arce et al. 



b|), and therefore we identify several molecular outflows in the 



L1641-N cluster using the ^^CQ ( J = 1 - 0) data 
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Molecular outflow activity in the L1641-N region was flrst found b y iFukui et a. 
at the position of the near-IR bright source IRAS 05338-0624 (see also 



' Fukui et a 


( 


1986) 


Fukui et al. 


1988 


)■ 



Furt her outflow surveys have been made through ra o. 



le.g., 

12 



Wilking et al. 



1990 



Stanke fc Williams 



20071 ). 



ecular line studies by seve ral authors 



Stanke fc Williams 



(120071 ) conducted 



CO (J = 2 — 1) mapping observations toward the L1641-N cluster using the IRAM 30 
m telesco pe and SMA inter: 



Recently, 



Davis et al 



'erometer, and identifled a number of molecular outflow lobes. 
((20091) carried out a near-infrared H2 survey of outflows in Orion 
A including the L1641-N region. To identify the molecular outflows, we flrst scrutinize 
the velocity channel maps of our ^^CO ( J = 1 — 0) (Figure E]) and the position-velocity 
diagrams to flnd localized blueshifted and redshifted emission. 

In Figure [m we present the ^^CO (J =1 — 0) integrated inte nsity contours overlaid on 



the 1.1 mm continuum image obtained by 



Shimajiri et al. 



(120111 ). The velocity intervals 



for the integration are ^ 5 km s ^ and 10 ^ 16 km s ^ for blueshifted and redshifted 



Stanke fc Williams 



components, respectively. 

from ^^CO (J = 2 — 1) whose names a re also designated iri 



coUimated redshifted lobe identifled by 



( I2OO7I ) labeled the outflow lobes identifled 
igur e im The most extended 



Stanke fc Williams 



( l2007l ). R-S, is not so prominent 



in our ^^CO (J = 1 — 0) map where this lo be is divided into s everal distinct knots. A 



number of knots labeled with SMZ 49 in 



Davis et al. 



(120091 ) are associated along this 



redshifted lobe. The blueshifted emission elongated along the northeast to southwest line 
is also seen north of the dust peak near the cluster center. This corresponds to the B-NE 



and B-N lobes in 



Stanke fc Williams 



(120071). In the ^^CO (J = 2 - 1) map, this lobe 



has a Y shape, whereas in our ^^CO (J = 1 — 0) map, the B-NE emission dominates 
over the faint B-N emis sion. This ^ ^ CO ( J =1 — 0) blueshifted lobe also corresponds 



Fukui et al. 



(119861 ). These are probably due to the fact that our 



to the one identifled by 
^^CO (J = 1 — 0) data tend to t r ace re latively low-velocity outflow components. In fact. 



according to 



Stanke fc Williams 



(|2007| ). at low velocities, the two lobes identifled by the 
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^^CO (J = 2 — 1) emission blend together with a broad patch of ^^CO emission. The 
most intense dust continuum emission is associated with the area between the redshifted 



(R-S) and blueshifted (B- NE and B-N 



and HCO+ ridge found by 



Fukui et aL 



obes , which is in good agreement with the HCN 



(119881 ) and the position of IRAS 05338-0624. The 
H2 flows labeled with SMZ 51 and SMZ 53 appear to follow the B-NE and B-N lobes, 



respectively (see Fig. 6 of IDavis et al. 



20091) 



Another highly-coUimated, redshifted lobe is see n to the west of the R-S lob e, which is 



(120071 ). Similarly 



in good agreement with the R-SW lobe in Figure 2 of lStanke fc Williams 
to the ^^CO (J = 2 — 1) map, this lobe appears not to h ave any obvious blueshifted 



counterlobe. The R-W, B-E, and B-SE2 lobes identified by 



Stanke fc Williams 



(120071 1 can 



be clearly recogni zed in our map. 



lobe identified by 



Fukui et al. 



mm . 



'he R-W lobe probably corresponds to the redshifted 



From our ^^CO map, the total outflow mass, momentum, and energy are estimated 
to about 13 Mq, 80 Mq km s~^, and 273 Mq km^ s~^, respectively, by assuming the 
local thermodynamical equilibrium (LTE) condition and optically-thin emission. Here, we 
adopt the same systemic velocity of 7.5 km s~^ f or all the outflow lobe s and the values 



are corrected for the inclination angle ^ = 57.3° (iBontemps et al. 



19961 ). The excitation 



temperature is adopted to 30 K, which is close to the peak brightness temperature of the 
^^CO ( J = 1 — 0) line profile at the L 1641-N cluster center. We also integrated the physical 
quantities in the velocity intervals of ~ 5 km s^^ and 10 ~ 16 km s~^ for the blueshifted 
and redshifted components, respectively. We note that the physical quantities estimated 
above are insensitive to the assumed excitation temperature. The estimated quantities 
increase only by a few xlO % over the range of Tex = 20 — 50 K. 

Our ^^CO map suggests that we identified the outflow lobes driven by about five 
embedded sources in the central part of L1641-N, the mean outflow momentum for a single 
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outflow may be estimated to be about 80/5 = I6M0 km s~^. If we assume tlie median 
stellar mass of 0.5 Mq, then this gives the outflow momentum per unit stellar mass of about 
32 km s~^, corresponding to the nondimensional parameter / of about 0.32, which gauges 
the strength of an outflow, where the wind vel ocity of 100 km s~^ is ado p ted. This value is 



consis tent with the fiducial ya.lues adopted by 



fl2006h . and 



Nakamura fc Li 



(I2OO7I I 



Matnzer fc McKeel (I2000l ) 



Li &: Nakamura 



5. Discussion 

5.1. Triggered Formation of Dust Filaments 

As shown in Section HJ some dust filaments are located near the cloud edge where the 
^^CO emission drops sharply, rather than along the ridge of the parent filamentary cloud 
(i.e., filaments A and E). The dust continuum emission associated with the filaments also 
tend to have very steep gradients in the envelopes of the filaments, in particular for the 
filament B. In addition, the dust filaments have two-components with different velocities. 
The velocity difference between the two components is significantly larger t han the typical 
local turbulent speed. According to the line-width-size relation obtained by Heyer Sz Brunt 
(12004J ). the typical turbulent velocity is estimated to be around 1 km s"\ about 3 times 
smaller than the velocity difference observed in the dust filaments. Thus, it is unlikely that 
the filaments were created by the dynamical compression due to the local turbulent flow. 
From these observational facts, we here propose that the dust filaments in this region were 
created by external compression, instead of the spontaneous gravitational contraction. 

On the larger scale of around 10 pc, there is a significant velocity gradient along the 
parent filamentary cloud. The cloud component with velocities smaller than about 6 km s~^ 
is dominant in the southern part of the parent filamentary cloud, whereas the component 
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with velocities larger than about 7 km is dominant in the northern part (see Figure |6] 
and Figure 2 of Bally et al. 1987). The interaction between the two- components is likely 
to have created the filaments A, B, C, and D in this region. Here, we refer to such a 
possible dynamical interaction as a cloud-cloud collision. Since the two components tend 
to overlap toward the dust filaments on the plane of the sky (see Figures [7|(b), [7](c), Mjo), 
and[H](c)), the collision may have occurred almost along the line of sight. The collision 
may also have triggered star formation in the L1641-N cluster as well as the formation of 
the massive filament C. The cluster formation triggered by a cloud-cloud collision hav e 



been recently discussed by several authors for other star- fo r ming regions (IXue fc Wu 



2008; 


Furukawa et al. 


2009 


Higuchi et al. 


2010 


)• 



Duarte-Cabral et al. 



2010 



2011 



Galvan-Madrid et al. 



2010 



Besides the possible cloud-cloud collision, other parts of this region appear to have a 
velocity structure that is in good agreement with an expanding motion created by a shell. 
A typical example is the filament E that has two velocity components inside the shell Al 
but appears to converge into a single velocity component at the outer part (see Figures [Tl^e) 
andlU^e)). Such an arc-like structure in the position- velocity diagrams can be interpreted 
as an expanding motion of the shell. The filament C, the most massive filament, also has 
such an arc-like structure in the position- velocity diagram in the western side (Figure [Sl^c)). 
This massive filament may have recently been compressed by an expanding shell after 
it was created by the cloud-cloud collision and the cluster formation was initiated. The 
compression due to the expanding shell may have accelerated the recent star formation in 
the L1641-N cluster. We don't think that the filaments A, B, C, and D were created by the 
expanding shells created by the stellar feedback from the protoclusters because their main 
axes tend to cross the shells and they appear more or less straight, rather than curved like 
the filament E. 
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Hence, the cloud-cloud collision and the expanding shells appear to have influenced the 
cloud structure and kinematics in this region. In Sections 15.21 and 15. 3^ we discuss these two 
dynamical events, the cloud-cloud collision and the expanding shell, in detail. 



5.2. Cloud-Cloud Collision 



According to our CO channel maps, the two-components with velocities of 4 — 6 km s~^ 
and 7 — 12 km s~^ appear to be interacting to form the dense filaments. The interpretation 
is consistent with the previous ^'^CO (J = 1 — 0) observations which suggested that two 
or more differe nt velocity component s are associa ted with the Orion A giant molecular 



cloud complex ( iBally fc Langer 



1987 



Bally 



19891 ). The velocity difference between the 



two-components is about 3 km s~^, corresponding to the Mach number of 10 under 
the assumption that the gas temperature is 20 K. The supersonic collision between the 
two-components would have increased the local density by a factor of 10^ if the shock is 
isothermal. Therefore, such a shock can create the dense filaments whose densities are as 
large as 10^~^ cm~'^ from preshock gas with 10^"'^ cm~^. In addition, the shock crossing 
time is estimated to be a few xlO^ — 10^ yr when the typical size of the filaments is adopted 
as a few pc. The estimated shoc k crossing time is s hort er than the lifetime of the L1641-N 



cluster (a few xlO^ yr, see, e.g.. 



Hodapp fc Deand (Il993[ ) and 



Galfalk fc Olofssonl (120081 )) 



by a factor of a few. Consequently, it is possible for such a collision to have triggered star 
formation in the dense filaments of this area. Since the shock crossing time is somewhat 
shorter than the lifetime of the cluster, the filament is expected to be in the postshock 
stage in which the observ ed two componeri t s may be passing over, thus leaving away from 



each other. According to 



Sakamoto et al. 



( 1l997l ). the blueshifted component is extended 



over the larger low-density area toward the eastern part of the cloud. This fact suggests 
that the collision with large clouds or flows (1—10 pc), instead of the collision with smaller 
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clouds, may have occ urred. The dynamical i nteraction with such external flow has been 
recently suggested by 



Shimajiri et al. 



( 120111 ) for the northern area of the Orion A giant 
molecular cloud complex. The diffuse blu eshifted component fou nd from our ^^CO data 



may be related to the flow pointed out by 



Shimajiri et al. 



(1201 ih . 



Recently, an interesting s cenario on the formation of t 



complex has been proposed by 



Hartmann &: Burkert 



le Orion A giant molecular cloud 



(120071 ) who performed two-dimensional 



numerical simulations of gravitational collapse of a thin gas sheet. According to their model, 
two large filaments are first created near the edge of the gas sheet by the gravitational effects 
of the cloud edge where the gravitational potential takes its local minimum (see also Larson 
1976 and Bastien 1983 for the effect of the cloud edge on the gravitational fragmentation). 
They demonstrated how an elongated rotating gas sheet with a density gradient along 
the major axis can gravitationally collapse to produce a structure qualitatively resembling 
the whole Orion A giant molecular cloud complex, having a fan-shaped structure at the 
southern part, ridges along the fan, and a narrow integral shaped filament at the northern 
part. In their model, our observed area, L1641-N, is located at the intersection between the 
fan and the narrow main filament, where the two large filaments collide with each other. 
Such a large scale interaction might be able to explain the origin of the dust filaments, 
although it remains unclear how the hypothetical parent gas sheet was created. 



5.3. Parsec-Scale CO Shells and Protocluster Winds 

As shown in Figure [9], our CO map revealed that several shells appear to shape the 
cloud structure in this region. The radii of the circles are about 10'— 20', corresponding to 
about 1—2 pc at a distance of 400 pc. Interestingly, the centers of the circles Al, A2, and 
A3 are close to the L1641-N cluster center, and those of the circles Bl, B2, and B3 to the V 
380 Ori center. Such large shells are unlikely to form by stellar feedback from a single young 
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star unless it is massive. Since the spectral types of the most massive cluster members are 
late B or early A for the two clusters, the shells may be difficult to create only by the stellar 
feedback of the most massive stars. Therefore, we suggest that these parsec-scale shells 
were created by multiple winds and/or outflows from cluster member YSOs. The reason 
why three or more shells are associated with each cluster remains unclear. One possibility 
is that the total momentum injection rate from cluster member YSOs is not constant 
but episodic because the star formation rate is likely to fluctuate or oscillate with time. 
Another possibility is that these circles represent dense parts of an expanding shell that 
is propagating into inhomogenius media Here, we call these YSO winds as "protocluster 
wind." In the next subsection, we consider the dynamical evolution of the shell driven by 
the protocluster wind, using a simple analytic model. 

Kinematical evidence that the protocluster wind is responsible for the cloud dynamical 
evolution comes from the filament E that is located at the western edge of the cloud and 
has an arc-like shape. The filament E has two velocity components inside the circle Al but 
appears to converge into a single velocity component at the outer part (see Figures [7](e) and 
m^e)), as mentioned in Section [STTl It is also distributed along the circle Al whose center 
coincides with the L1641-N cluster, suggesting that it may have been created by the shell 
driven by the protocluster wind from the L1641-N cluster. 

Another example of the parsec-scale expanding shell is the one whose center is close 
to V 380 Ori. This s hell is fully covered by the ^"^CO (J = 1 — 0) map obtained by 



Bally fc Langer 



( ll987l ). although the southern part of the shell is not covered by our ^^CO 
data. We show in the ^^CO integrated intensity map presented in Figure [12] the circles that 
appear to fit the shells. The ^^CO position-velocity diagrams of the shells toward the two 
positions (Figures [^d) and[8t^f)) clearly show the existence of the two different velocity 
components, consistent with the expanding motion. The expanding motion of the shell is 



in good agreement with iHeyer et al. I ( 119921 ) who found a similar expanding motion of the 
shelL 



5.4. Expanding Motions Driven by Protocluster Winds 



In the following, we discuss how the protocluster wind evolves in the parent molecular 
cloud using a simple analytic model. Here, we consider the motion of an expanding shell 
driven by protostellar winds from cluster member YSOs in a uniform media with density 
Pq. For simplicity, we don't take into account the effects of protostellar outflows in the 
following calculations, but discuss the effect of the protostellar outflow feedback in the next 
subsection. 

If a wind from a protostar with a radius i?* is driven by a ram pressure PwV^, then the 
motion of an expanding shell is described as 



where p„ is the wind density, Kv is the wind velocity, R is the radius of the shell, and N is 
the number of YSOs. For simplicity, all the stars are assumed to inject the wind momentum 
simultaneously at the same constant rate. From the above equation, the time evolution of 
the shell radius is given by 




(1) 



R 




1/2 




M 



) 



1/4 



) 



1/4 



~ 2 pc 



w 



-1 



200 km s 



1 




(2) 



where M is the mass loss rate from a single YSO and given by AirRlpy^Vy^. The expanding 
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dR 
lit 



velocity of the shell, K, is evaluated as 
Vs = 

^ 1 km s~^ 
Po 



-) 

100 J 



1/4 



M 



1/4 



lO-Woyr- 



200 km s" 



1/4 



-1/4 



t 



-1/2 



(3) 



103 cm-3/ \^106 yr 

The typical mass loss rate and wind velocity from YSOs are somewhat uncertain. Some 
observations suggest that the mass loss rate and wind velocity from a low-mass YSO are o f 



the order of 10 "^Mc^ yr ^ and a few xlO^ km s ^, respectively (e.g. 



Wilkin fc Stahler 



Norman fc Silk 



1980; 



1998h . Here, we adopted M = 10''^ Mq yr'^ and V„ = 200 km s'^ as 



their representative values. For the L1641-N and V 380 Ori clusters, the number of cluster 
member YSOs is of the order of 100. Adopting the density of 10^ cm"^ for the ambient 
gas and the cluster lifetimes of 10^ yr, the radius of the expanding shell and the expansion 
velocity are evaluated to be about 2 pc and 1 km s~^, respectively, comparable to those of 
the shells identified from our ^^CO data. Therefore, the protocluster winds from these two 
clusters are likely to have enough energies to produce the parsec-scale shells within a few 
Myr. 



5.5. Molecular Outflows in the L1641N cluster 



In the previous subsection, we omitted the contribution of the outflow feedback on 
the dynamical evolution of the expanding shell. However, recent numerical simulations 
of cluster formation have suggested that protostellar outflow feedback plays an important 
role in regulating star formation in cluster-formi ng clumps because they iniect substantia l 



amount of kinetic energy into the surround ings (ILi fc Nakamura 



2007 



2006 



Nakamura fc Li 



Carroll et al. 



2009 



Wang et al.l 120101) . In addition, the propagation directions of 
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the collimated outflows in the clustered environment are not so aligned preferentially 
i n the global magnetic field direction even in the presence of a strong magnetic field 



( iNakamura fc Li 



2011 



). Therefore, the total outflow momentum is expected to be injected 
isotropically on average. Here, we attempt to assess how the outflow feedback of the 
L1641-N cluster influences the surrounding gas and contributes to the dynamical evolution 



of the expanding shell. 



5.5.1. Dynamical State of the LI64I-N Cluster Forming Clump 



Recently, INakamura et al.l (l2011bl ) investigated how molecular outflow feedback 
influences the dynamical state of a nearby cluster-forming clump. Serpens South, by 
applying the virial analysis. They found that the Serpens South clump is close to virial 
equilibrium and the total kinetic energy injected by the current outflow activity is less than 
the clump gravitational energy, concluding that the current outflow acti vity is not enoug h 



to dest roy the whole cluster-forming clump. In the following, we follow 



Nakamura et al. 



(l2011bl ) and clarify how the current outflow activity in the L1641-N cluster impacts the 



whole clump. 

The virial equation for a spherical clump is given by 

where the terms, /, f/, and W , denote the moment of inertia, internal kinetic energy, and 
gravitational energy, respectively. Here, we neglect the surface pressure term. A clump is in 
virial equilibrium when the right-hand-side of Equation (jll) is zero, i.e., 2U + = 0. The 
kinetic energy and gravitational energy terms are expressed, respectively, as 

U^^-M^ (5) 
161n2 ^ ^ 
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and 



W = -a 



R 



1 



(6) 



where M is the clump mass, AV is the ID FWHM velocity width, G is the gravitational 
constant, R is the radius of the clump, the values $ and are, respectively, the magnetic 
flux penetrating the clump and the critical magnetic flux above which the magnetic field 
can support the clump against the self-gravity. Although according to recent observations 
of magnetic fields associated with the nearby p arsec-scale cluster-forming clumps, the valu es 



of ^/^rr are estimated t o be around 0.5 (e.g 



2011 



Kwon et al. 



Falgarone et al. 



2008 



Sugitani et al 



2010 



20111 ). we here neglect the effect of the magnetic support for simplicity 



i.e., $ = 0). The dimensionless parameter, a, measures the effects of a nonuniform and/or 



nonspherical mass distribution ( iBertoldi &: McKee 



19921 ). and is of the order unity. For a 
uniform sphere and a centrally-condensed sphere with p oc r~^, a = 3/5 and 1, respectively. 
Here, we adopt a = 1 because the cluster-forming clump tends to be centrally condensed. 



According to 



Stanke &: Williams 



(120071 ). the radius and mass of the cluster-forming 
clump associated with the L1641-N cluster are roughly estimated to be 0.5 pc and 150 Mq, 
respectively. Here, only the part involved in directly building the L1641-N cluster, i.e., a 
head of the clump, is defined as the cluster-forming clump, although the parent clump has 
a long tail extending southeast of the L1641-N cluster, as shown in the 1.1 mm map. The 
ID FWHM velocity width is estimated at about 2 km s~^ from the ^^CO data. Using these 
values, the kinetic energy and gravitational energy terms are estimated to be 2U 325M0 
km^ s~^ and W ~ — 194Mf7i krn ^ s"^, respectively. These values are comparable to the 



estimate by 



Stanke fc Williams 



(120071 ) who conducted a similar analysis using their ^^CO 
(J = 2 — 1) data. Thus, the clump is likely to be close to virial equilibrium or somewhat 
gravitationally unbound. On the other hand, the total kinetic energy due to the outflows is 
evaluated at E^^t ~ 273 Mq km^ s~^. This is comparable to the clump gravitational energy. 
Although the substantial amount of the outflow kinetic energy appears to escape from the 
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clump, the energy input due to the current outflow activity seems to influence the clump 
dynamics significantly. 

5.5.2. Effects of the Molecular Outflows on the Expanding Shells 

As shown in Figure [111 the identified outflow lobes often extend beyond the dense 
clump traced by the dust continuum emission, implying that the outflow momentum 
escaped from the clump is likely to influence the density and velocity structure in the 
surroundings. If the escaped momentum fllux is about a half the total injected outflow 
momentum flux, it is estimated to be about 1Q~^Mq km s^^ yr^^, assuming the outflow 
dynamical time of a few xlO^ yr. On the other hand, the total momentum flux injected 
from the protostellar winds is estimated to be NMV^ ~ 100 x lO"'' x 200 = 2 x lO^^M© 
km s^^ yr^^, where we adopted the values given in Equations ([2]) and ([3]). The estimated 
wind momentum flux is comparable to that of the protostellar outflow feedback. Therefore, 
we conclude that both the protostellar outflow and wind feedback can contribute to the 
dynamical evolution of the expanding shell created by the protocluster wind from the 
L1641-N cluster. However, the weak dependence on the injected momentum in Equations 
([2]) and ([2D indicates that even if we take into account both the protostellar outflows and 
winds, the estimated radius and expanding velocity of the shell do not change significantly. 

6. Summary 

We carried out the ^^CO {J = 1 — 0) mapping observations toward the 48' x 48' area 
including the L1641-N cluster in the Orion A giant molecular cloud complex, using the 
Nobeyama 45 m telescope. The main results are summarized as follows. 

1. From comparison between the ^^CO (J = 1 — 0) map and the 1.1 mm continuum 
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image, we found that several dust filaments are located near the cloud edge traced by the 
i^CO ( J = 1 - 0) emission. 

2. The dust filaments have two-components with different velocities. The velocity 
difference between the two-components is about 3 km s"-*^, which is significantly larger than 
the typical local turbulent speed of 1 km s~^. Therefore, we suggest that the dust filaments 
were created by the dynamical compression triggered externally, i.e., a cloud-cloud collision, 
instead of the spontaneous gravitational contraction. 

3. The ^^CO ( J = 1 - 0) and ^^CO ( J = 1 - 0) velocity channel maps suggest that the 
blueshifted (Vlsr < 6 km s~^) and redshifted (Vlsr ^ 7 km s""^) components are interacting 
with each other. Since the two components appear to overlap toward the dust filaments on 
the plane of the sky, the collision between the two components may have occurred almost 
along the line of sight. A good example of such a cloud-cloud collision is the most massive 
dust filament, associated with the L1641-N cluster, which has the two-components with 
different velocities at the long tail stretching from the head of the dust filament. We suggest 
that the formation of the L1641-N cluster may have been triggered by such a coUision. 
Since the shock crossing time is somewhat shorter than the lifetime of the L1641-N cluster, 
the massive filament associated with the cluster may be in the postshock stage in which the 
two components with the different velocities are passing over, thus leaving away from each 
other. 

4. We found several parsec-scale shells in the ^^CO {J — 1 — 0) data cube. Some of the 
shells appear to be spatially well-ordered and homocentric. The centers of the shells are 
close to either the L 1641-N or V 380 Ori cluster centers, implying that the star formation 
activity in the clusters may be responsible for the formation and evolution of the shells. In 
particular, the shell surrounding V 380 Ori is prominent in the ^^CO map. 

5. The molecular gas distribution and kinematical structure of this region led us 
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to the following scenario. On the large scale of at least about 1 — 10 pc, a cloud-cloud 
collision may have occurred almost along the line of sight in this region, contributing to 
the formation of several dense filaments (the filaments A through D). The cloud-cloud 
collision have triggered the formation of the L1641-N cluster. Multiple protostellar winds 
and outflows from the cluster member YSOs created large expanding bubbles that can be 
recognized in the ^^CO and ^'^CO maps. Here, we call these YSO winds as "protocluster 
winds." The dust filament located at the western cloud edge (filament E) appears to 
curve along the shell whose center is the L1641-N cluster, and have two different velocity 
components at the inner part of the shell. The two components appear to converge into 
a single velocity component at the outer part of the shell. It is consistent with an idea 
that the filament E is a part of the expanding shell created by the protocluster wind 
from the L1641-N cluster. The shell surrounding V 380 Ori also has two different velocity 
components both in the ^^CO and ^'^CO maps, reaching the filament C. It presumably hit 
the filament C, influencing the recent star formation in the L1641-N cluster. Both the 
cloud-cloud collision and the protocluster winds are likely to have created the complicated 
cloud morphology and kinematics in this region. 

6. Toward the L1641-N cluster, we identified a number of the outflow lobes using 
the ^^CO (J = 1 — 0) data. The i dentifi ed lobes are reasonably in good agreement with 



Stanke &: Williams 



( 120071 ) who identified the outflow lobes in this region 



the results of 

with the ^^CO (J = 2 — 1) emission, using the IRAM 30 m telescope and SMA. The 
total outflow energy in the L1641-N cluster is comparable to the gravitational energy of 
the cluster-forming clump. This may suggest that the outflow feedback influences the 
dynamical evolution of the clump signiflcantly. Assuming the median stellar mass of 0.5Mq, 
the mean outflow momentum per unit stellar mass is estimated to be about 32 km s~^, 
under the assumption of optically thin gas. This mean outflow momentum corresponds to 
the nondimensional outflow parameter of / = 0.32, which gauges the strength of an outflow. 
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Li fc Nakamura 
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Nakamura &: Li 
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Fig. 1.- 



12 



CO {J = 1 — 0) total integrated intensity map in the velocity range from 



fLSR = 0.0 to 20.0 km toward the 48' x 48' area including the L1641-N cluster. The 
effective spatial resolution of the map is 21". The ^^CO emission peak coincides reasonably 
well with the position of the L1641-N cluster, (b) Same as panel (a), but the 1.1 mm 
continuum data are overlaid on the ^^CO (J = 1 — 0) total integrated intensity map. The 
contours are drawn at 0.06, 0.15, 0.3, 0.6, and 1. 2 Jy/beam. The l a rm s noise level is 9 



mJy/ 



RebuU et al. 



(120061) and 



Carpenter et al. 



3eam. The positions of YSOs identified by 
( 1200 ll ) are overlaid on the images by the circles and crosses, respectively. Note that the 
plotted YSOs are less embedded sources, and therefore younger embedded YSOs are not 
shown. Several dust filaments are labeled with alphabet A through E. The positions of 



L1641-N and V 380 Ori are indicated by the white arrows. 
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Fig. 2. — (a) ^'^CO {J = 1—0) total integrated intensity map toward the same area presented 
in Figure [T]in the ve l ocity range from wlsr = 1-0 to 14.0 km s^^. The ^^CO data were obtained 



by 



Bally fc Langerl (119871 ) with the 7 m AT&T Bell laboratories telescope. The strongest 



^■'CO emission is associated with the L1641-N cluster, (b) Same as panel (a) but the 1.1 
mm continuum data are overlaid with the contours on the ^^CO {J = 1 — 0) total integrated 
intensity map. The contours are drawn at 0.06, 0.15, 0.3, 0.6, and 1.2 Jy/beam. Several 
dust filaments are labeled with alphabet A through E. 
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Fig. 3.- 



12 



CO {J = 1 — 0) peak intensity map. (b) Same as panel (a) but the 1.1 



mm continuum contours are overlaid. The contours are drawn at 0.06, 0.15, 0.3, 0.6, and 
.2 Jy/beam. The crosses and circles defnotes the positions of YSOs identified by Spitzer 



( iRebuU et al. 



20061), respectively. Note that these YSOs are less embedded sources. Several 



dust filaments are labeled with alphabet A through E. 
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Fig. 4. — H^'^CO^ ( J = 1 — 0) velocity integrated intensity contours overlaid on the 1.1 mm 
dust continuum image taken by AzTEC on ASTE. The velocity range is from wlsr = 4.0 to 
14.0 km s~^. The contours start from 0.42 K km s~^ at intervals of 0.2 K km s"^. The solid 
lines indicate the observation box for the H^^CO^ {J = 1 — 0) emission. The data were taken 



from the Nobeyama 45 m archival data (http://www.nro.nao.ac.jp/ ). Several dust filaments 
are detected in the H^^CO^ emission, indicating that they contain dense gas with 10^ cm~^. 
Several dust filaments are labeled by alphabet A through E. 
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Fig. 5. — ^^CO ( J = 1 — 0) velocity channel maps with velocity width of 1.0 km s ^. 
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Fig. 6. — Two-color image toward the same area presented in Figure [T] with redshifted ^^CO 
(J = 1 — 0) integrated intensity in red (9.5 km s"^ ~ 14.5 km s~^) and blueshifted ^^CO 
(J = 1 — 0) integrated intensity in blue (3.5 km s^^ ~ 6.5 km s~^). The white contours 
indicate the 1.1 mm dust continuum emission at levels of 0.06, 0.15, 0.3, 0.6, and 1.2 Jy/beam. 



- 39 - 




Fig. 7. — (a) ^^CO {J = 1—0) peak intensity map showing the positions of the PV diagrams 
presented in panels (b) through (e). The 1.1 mm continuum map is overlaid on the image 
with the black contours whose levels are the same as those of Figure [TJ (b) Position- velocity 
diagram along the line (1). The abscissa denotes the offset measured from the position of 
the white cross indicated in panel (a). The plus and minus values are for the western and 
eastern sides of the cross, (c) Same as panel (b) but for the line (2). (d) Same as panel (b) 
but for the line (3). (e) Same as panel (b) but for the line (4). 
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13 



CO (J 



0) velocity integrated intensity map showing the positions of 



the PV diagrams. The velocity integration range is the same as that of El The lines (b) 
through (e) are the same as those of Figure [71 The 1 . 1 mm continuum map is overlaid on 
the image with the black contours whose levels are the same as those of Figure [H Note 
that the southern part (be low Dec. ~ —6 : 30) was not covered by the 1.1 mm observations 



see 



Shimajiri et al. 



2OIII ) (b) Position- velocity diagram along the line (1). The velocity 



resolution of the data is 0.2 km s~^. (c) Same as panel (b) but for the line (2). (d) Same as 
panel (b) but for the line (3). (e) Same as panel (b) but for the line (4). (f) Same as panel 
(b) but for the line (5). 



-41 - 




Fig. 9. — Same as Figure [3] but the positions of parsec-scale shells are indicated with the 
dashed lines. 



Fig. 10. — Three-dimensional representation of the antenna temperature (T^) of the CO 
(J = 1 — 0) emission in the RA-DEC-l^sR space. The green color roughly represents the 
parts with ~ 15 K. Several arc-like structures can be recognized in the RA-DEC-Vlsr 
space. 
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Fig. 11.— CO (J = 1-0) 

1.1 mm image obtained by 



integrated intens i ty con tours toward the L1641-N cluster on the 



limaiiri et ah 



as that shown in Figure 2 of 



(2011 



Stanke fc Wilhams 



The displayed area is almost the same 



( 120071 ) who mapped the same area in CO 



(J = 2 — 1). The outflow lobes seen in the CO ( J = 1 — 0) data are also labeled following 



Stanke fc Williams 



( 120071 ). The blue contours represent blueshifted CO (J = 1—0) intensity 
integrated from 0.25 km s^^ to 5.25 km s~^, starting from 5 K km s~^ at intervals of 0.75 
K km s~^. The red contours represent blueshifted CO ( J = 1 — 0) intensity integrated from 
9.75 km s~^ to 15.75 km s~^, starting from 9 K km s~^ at intervals of 0.75 K km s~^. The 
gray scale represents the 1.1 mm image in units of Jy beam~ ^. The positions of the classical 



and weak-line T-Tauri stars observed by 
crosses, respectively. 



Fang et al. 



(120091) are shown in black and green 
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Fig. 12. — (a) ^^CQ tota 
km (IBally fc Langeii 



■ integ rated intensity map in the range from Vlsr = 1 km s ^ to 14 
19871). (b) 13 CO peak intensity map in the same area presented in 



panel (a). For both the panels, the positions of the circles Al and Bl presented in Figure [9] 
are indicated with the dashed lines. 



